Search for universality in one-dimensional ballistic annihilation kinetics. 
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We study the kinetics of ballistic annihilation for a one-dimensional ideal gas with continuous 
velocity distribution. A dynamical scaling theory for the long time behavior of the system is derived. 
Its validity is supported by extensive numerical simulations for several velocity distributions. This 
leads us to the conjecture that all the continuous velocity distributions 4i{v) which are symmetric, 
regular and such that (j){0) ^ are attracted in the long time regime towards the same Gaussian 
distribution and thus belong to the same universality class. Moreover, it is found that the particle 
density decays as n{t) ~ t'", with a ~ 0.785 ± 0.005. 

PACS numbers; 82-20. Mj, 05.20Dd 



I. INTRODUCTION 

Ballistically-controUed reactions provide simple exam- 
ples of non-equilibrium systems with complex kinetics 
and have recently attracted a lot of interest They 
consist of an assembly of particles moving freely between 
collisions with given velocities. When two particles meet, 
they instantaneously annihilate each other and disappear 
from the system. In one dimension, it is enough to con- 
sider point particles and we shall restrict ourselves to 
this case here. The system with only two possible ve- 
locities +c or — c has been studied in a pioneering work 
by Elskens and Frisch Using combinatorial analy- 
sis, they showed that, in the long time limit, the den- 
sity of particles was decreasing according to a power law 
in the case of a symmetric initial velocity distribu- 
tion; Krug and Spohn Q obtained independently similar 
results. Later on, Redner et al. |^,^ have studied the 
case of more general velocity distributions. Based on nu- 
merical simulations and mean-field type arguments they 
showed that the exponent characterizing the power law 
decay of the particle density could depend on the veloc- 
ity distribution. The case of a general distribution has 
been recently approached analytically by Piasecki and 
Droz et al. ij^ji^l- It was shown that the annihilation dy- 
namics reduced exactly to a single closed equation for 
the two-particle conditional probability Q. A method 
which permits to solve this evolution equation for dis- 
crete velocity distributions was developed in Q and ex- 
plicitly applied to the case of a symmetric three- velocity 
distribution. It turns out that on one hand, different dy- 
namical behaviors can occur depending on the relative 
probability weights given to the three velocities and, on 
the other hand, that the fluctuations are playing a very 
important role, invalidating the predictions of mean-field 
or Boltzmann like approaches. 

The case of continuous velocity distributions was con- 
sidered first by Ben-Naim and co-workers 0]. Within a 



Boltzmann approximation, they showed that the parti- 
cle density was decreasing according to a power law 
Moreover, for initial velocity distributions with a power- 
law dependence near the origin and a cut-off at vq, i.e. 



e{v^ - \v\) (-1 < ^ < 0), they found that the 



exponent a was non-universal and was depending contin- 
uously on the value of /x. This dependence was confirmed 
by Monte-Carlo simulations. Note that the particular 
case /i = corresponds to a uniform distribution. In this 
case, and in one dimension, numerical integration of the 
Boltzmann equation led to a = 0.77, while Monte-Carlo 
simulations gave a = 0.76. These values are very dif- 
ferent from the ones obtained for the discrete two- and 
three-velocity distributions. 

The distributions considered by Redner et al. form a 
very special subset of the continuous velocity distribu- 
tions. So, it would not be justified to generalize their 
conclusions. A deeper understanding of the problem re- 
quires to study a wider class of distributions. This was 
our motivation for this work. Several questions could be 
asked. The main issue analyzed here, concerns the ex- 
istence of universality classes in this problem. Namely, 
does it exist a class of continuous velocity distributions 
for which all its members possess the same exponent a, 
and what are the characteristics of this class? As we 
have seen above, to approach this question it is not suffi- 
cient to consider a mean-field like approximation, but it 
is necessary to build up a theory taking into account the 
fiuctuations. 

This paper is organized as follows. In Section ^ we 
define the model and present the main ideas and conclu- 
sions obtained within the analytical approach developed 
by Piasecki and Droz et al. Based on those re- 

sults a dynamical scaling theory is developed in Section 
III , whose validity is confirmed by extensive numerical 
simulations in Section [V . Three different continuous ve- 



locity distributions (Gaussian, uniform and Lorentzian) 
are studied. We found that they are all attracted dur- 
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ing the time evolution towards a Gaussian distribution 
for which the decay exponent is a = 0.785 ± 0.005. It 
is conjectured that all the continuous velocity distribu- 
tions regular near the origin and having a finite nonzero 
value for v — are attracted by the same Gaussian dis- 
tribution and thus belong to the same universality class. 
Moreover, it is shown why the velocity distributions with 
a power-law dependence with negative exponent fi near 
the origin are not attracted by the Gaussian distribution 
and can lead to non-universal exponents. Final remarks 
are made in Section ^ and the numerical algorithm used 
in the simulations is explained in Appendix. 



II. THE MODEL 

We assume that initially the particles are uniformly 
distributed in space, according to the Poisson law, with- 
out any correlations between their velocities. Note that 
other distributions than Poisson could be considered as 



long as one is dealing with a renewal process (see 1 10 for 
the definition of a renewal process). The fundamental 
role in the analysis made in is played by the distri- 
bution of nearest neighbors. Suppose that at time t there 
is a particle at point xi in the fluid, moving with velocity 
vi. We denote by 



/J.{x2,V2\xi,Vi;t) 



(1) 



the conditional probability density for finding its right 
nearest neighbor at distance 3:21 — X2 — xi > 0, with 
velocity V2- The density (0) satisfies the normalization 
condition 

(i2^(2|l;i) = l (2) 

where a convenient short hand notation 

j = {xj,Vj), dj = dxjdvj, j = 1, 2, . . . 

has been used. As a rigorous consequence of the dynam- 
ics of ballistic annihilation , this quantity obeys the 
equation: 



at 0x1 0x2 J 

di [C(l, 3)At(3|l; t) - C(2, 3V(3|2; t)]^l[2\l■t) 

d3 [ d4C'(3,4)/i(3|l;t)/i(4|3;t)/i(2|4;t). (3) 



where C(l, 3) is the binary collision operator defined as: 

C{l,2)^vi20{vi2)Six2i-0+) (4) 

9{x) is the usual Heaviside function and V12 stands for 
vi ^ V2- If each particle has initially the same con- 
tinuous probability density (t>(v) to move with veloc- 
ity V, the state of the system is translationally invari- 
ant. As a consequence, /i(2|l;i) depends in the posi- 
tion space only on the distance X21 [we shall then write 
fi{x2,V2\xi,Vi;t) — /i(x2i, W2IO, Wi; <)]. A particular role 
is played by the value of density at contact 



fi{0^ ,V2\0,vi;t) = lim fi{x,V2\0,vi;t) 



(5) 



The notation 0+ stresses the fact that the distance be- 
tween the particles which are to collide approaches zero 
through positive values. The higher order conditional 
distributions /is(2, 3, . . . , s|l; t), s = 3, 4, . . . factorize for 
all times i > as 

fis{2,3,...,s\l;t) = l[fi{j\j-l;t), s = 3,4,... (6) 

providing that this condition is fulfilled at time t = 0. As 
a consequence, /i(2|l;t) has the extraordinary property 
to obey exactly the closed nonlinear integro-differential 
equation (||). 

As we shall see later, it will be more convenient to con- 
sider iy(x2,V2\xi^vi;t) = i'{2\l;t) the (conditional) den- 
sity of particles with velocity V2, which are at a distance 
X2 — xi of a given particle of velocity vi. With Eq. ^ it 
follows that iy{2\l;t) can be expressed as: 



i/(2|l;t) = ^(2|l;i) + ^ / d3 ■ ■ ■ f ds n{2\s;t)fi{s\s - l;t) . . . fi{3\l;t). 

s=3 •' 



(7) 



We can then write a new closed equation governing the dynamics of v which, for a translationary invariant system, 
takes the form: 



^ + "21^ + V2lO[Vi2)d(X 



iy{x,V2\0,vi;t) 

= I dv3{6{vi3)vi3v{0-^,V3\0,vi;t)[iy{x,V2\0,V3;t) ~ v{x,V2\^,v\\t)\ 



+ Q{v2z)v2-iv{Q'^ , walO, V2;t)v{x, W2 10, wl; t) 
+ Q{v-i2)vz2v{Q^ , W3IO, vi\t)v[x, t)2|0, i;3; t)}. 



(8) 
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A complete description of the dynamic of the system 
is obtained once the equation of motion for f{v;t), the 
density of particles with a velocity v, is given. It is read- 
ily obtained from the dynamics of ballistic annihilation. 
One finds: 



d_ 

dt 



f{vi\t) 



dv2 0{vi2)vi2v{0'^, V2\0, vi;t)f{vi;t) 

dV2 0(W21)W21I^(O+, Vi\0, V2]t)f{v2]t). (9) 



For a symmetric initial condition [i.e. (j){v) — 0(— w)], 
this last equation takes a simplified form 

— = / dv2 [6{vi2)vi2v{Q+ ,V2%vi-t) 



0(v2iKiKO+,-«2|O,~«i;i)], 



(10) 



where we used that for a symmetric distribution 

f{vi\t)v{x,V2\Q,Vi;t) = f{~V2;t)v{x,-Vi\Q,-V2\t). 

To get a closed analytical form for the particle density 
f{v;t), we thus have to solve Eq. (^), or at least to ob- 
tain the density at contact 1/(0^ ,V2\0,vi;t). This is a 
very difficult problem and despite numerous attempts, 
we have not yet been able to complete this program. Ac- 
cordingly, we choose a less ambitious approach. In the 
next section, we shall propose a scaling theory for f{v;t) 
and 2^(0"*", f2|0, ui; whose validity will be confirmed by 
numerical simulations. 

III. SCALING THEORY 

We are mainly interested in the behavior of the sys- 
tem in the long time limit, where the time dependence of 
the physical quantities of interest is given by power laws. 
Thus, it is natural to try to develop for them a dynamical 
scaling theory. The validity of this scaling theory will be 
attested a posteriori by numerical simulations. 

For simplicity, we shall restrict our analysis to sym- 
metric distributions. For dimensional reasons, a time 
dependent characteristic velocity should be introduced. 
For a symmetric distribution, the average velocity van- 
ishes for all times. However, the mean-square velocity, 
denoted M{t), does not vanish. Thus we shall assume 

1/2 

that [M(t)] ' is the time dependent characteristic ve- 
locity which will enter into the scaling theory. 

Hence, we postulate the two following scaling laws: 



[M{t)Y 



iy{0+,V2\0,vi;t) 



n{t) 



-1/2 



(11) 



X <i>{0+,V2[M{t)r'^^\0,v,[M{t)]-'^^y (12) 



and 



n{t)^ jdvf{v-t) 
Jdvvy{v;t) 



Jdvfiv;t) 



(13) 



(14) 



It is useful to introduce the intermediate quantity [see 
Eq. d)]: 



Ic{vi;t) = J dv2[9ivi2)vi2iy{0'^ ,V2\0,vi;t) 

+ 0{v2i)v2M0+ ,-V2\0,-vi;t)]. (15) 
Using the scaling forms (pi]) and (p2|), one can write: 

h{vi;t)^n{t)[M{t)f^ 



dU2 6'(ui2)wi2<l'(0+, M2IO, Ml) 

+ j du2 6I(m2i)m21$(0+, -U2IO, -Ui) 

= n{t)[M{t)f'^Uui), (16) 

where Ui = Vi[M{t)] , {i = 1,2). The time deriva- 
tive of the density and the mean square velocity can be 
written as: 



n{t) ^ Jdvl,{v;t)f{v;t) 
n{t) Jdvf{v;t) 

= -nit)[Mit)] 



1/2 J duJc{u)^iu) 



Jdu^^ii 



(17) 



M(i) ^ n{t) _ Jdvv^Ic{v;t)fiv;t) 



M{t) n{t) 



J dv u^/(f ; t) 



/^^^«-r/^Ml/2 /^" uVc(m)^(m) 

-n{t)[M{t)] jr- -— — . (18 

J du M^ w(m) 



We note that Jdw Jc(u)^('")//rfw ^(u) and Jduu^Jc{u) 
xv['(it)/J duu^'i/{u) are two constants called A and B re- 
spectively. In the long time limit, the solution of Eqs 
(0) and (|l|) is: 



n(t) - t-°' 

M{t) - r'^i^ 



with 



a = 2A/{A + B) 
P={B-A)I{A + B). 

For all values of A and B, the scaling law: 

a + /3= 1 



(19) 
(20) 



(21) 
(22) 



(23) 



holds. Note that this scaling law has already been ob- 
tained by Ben-Naim et al. from heuristic arguments. 

Using this scaling law, it is possible to justify a poste- 
riori our scaling postulate for v. Indeed, let us introduce 
the dimcnsionless correlation function g{x,V2,vi;t) sat- 
isfying 



where n{t) is the particle density at time t: 



iy{x,V2\0,vi;t) = g{x,V2,vi;t)f{v2;t). 



(24) 
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X2 





X2 - V^t 



t = 



b) 




X2 - V2t) 



bt 



t = 



FIG. 1. Illustration of the similarity relation satisfied 
by the correlation function g{x2i,V2,vi;t). a) A typical 
two-particle configuration is shown at time t. Initially the 
positions of the particles are xi — vit and X2 — V2t. They 
move respectively with velocity vi and V2. b Same configu- 
ration in which time and distances have been rescaled by a 
factor b (the velocities remaining unchanged). The trapezoid 
defined by the points b{xi — vit), b{x2 — V2t) at t = 0, and 
bx2, bxi at bt in Fig. b) is similar to the one defined by the 
points xi — vit, X2 — V2t at t = 0, and X2, x\ at t in Fig a). 

As shown on Fig. |l|, by rescaling both length and time 
by a factor 6, g should satisfy the following similarity 
relation: 

g(x,V2,vx;t) = g{bx,V2,vi;bt). 

Moreover, rescaling the velocities by a factor c corre- 
sponds to a dilatation of the time scale by the same fac- 
tor: 

g{x,cv2,cvi\t) = g{x,V2,vi]ct). 

Combining these two relations and putting h = n(t) and 
c = [M(t)]^^^, one immediately finds 

/ —1/2 

g{x,V2,vi;t) = g\n{t)x,V2[M{t)] , 

v,[M{t)]-'/^;n{t)[M{t)f't). (25) 

For sufficiently long time, the scaling law (^) implies 

that n(i)[M(t)]^/^ - t-\ Using Eq. (||) and the scal- 
ing form ( pi] ) for /, one recovers then our scaling postu- 
late (p^. This confirms the self-consistency of our the- 
ory. Furthermore by introducing the scaling postulates 
(|ll|,|2|) in (0) one obtains: 



3A- B 



B~A 



u*'(u) 



hence an expression relating '^{u) to Jc{u): 
3A- B 



B-A *'(u) 



Jc(u). 



(26) 



(27) 



Let us now suppose that, initially, /(u; 0) 



\v\), with — 1 < /i < 0. In this case Ben-Naim et al. |^ 
argue that the long time behavior of n{t) is nonuniversal 
(/i-dependent) and that ^(u) retains the same power-law 



dependence than /; in particular it diverges like 



for 



M — + 0. What could we say about this case within our 
scaling approach? Letting w ^ in Eq. (p7|), we find 



3A- B B-A 



or in terms of the exponent a: 



2a -I 
1 - a 



2Jc(0) 
B-A' 



(28) 



(29) 



Jc(0) might in principle depend on /i. However, it is very 
unlikely that its dependence will be such as to exactly 
compensate the linear term in /z, hence the nonuniversal- 
ity of the exponent a. On the contrary, if \E'(m) satisfies 



lim u—^ = 



(30) 



[which is the case if ^(u) and ^''(u) are continuous func- 
tions at u = and if ^'(0) 7^ 0], wc find that 



2a - 1 _ 2Jc(0) 
1-a ^ B-A' 



(31) 



Hence a will be universal only if the right hand side does 
not depend on the details of the initial distribution. We 
are unfortunately unable to prove analytically this inde- 
pendence on the initial conditions. However, in the next 
section we shall show that this property is supported by 
precise numerical simulations for several initial velocity 
distributions. Furthermore we shall find on the same 
footing that 5" (it) is an universal Gaussian function. Note 
that for ^'(u) to be universal, one should have according 
to Eqs (11) and (|l|) that Jc{u)/{B - A) is universal. 



IV. NUMERICAL SIMULATIONS 



To test the validity of our scaling forms (11) and (|lj) , 
we have performed numerical simulations in one dimen- 
sion. The method we used is an exact synchronous time 
evolution whose algorithm is detailed in Appendix. We 
considered the following three different initial continuous 
velocity distributions: 



i) A Gaussian distribution: 



1 



27rwo 



exp 



(32) 



with vq = [M(0)]^/^ = 1. 
ii) A uniform distribution with a cut-off: 

(l^{v)^:^e{vo-\v\), 
with vo — V3, such that M{0) ~ 1; 



(33) 



4 



iii) A Lorentzian distribution 



with characteristic velocity vq — I. 

For each simulation, we started with 2^^ particles on a 
periodic chain of length 2 . Note that the initial char- 
acteristic collision time Tc = [n(0)wo]~^ is of the order 
10~^. We choose to distribute the particles uniformly on 
the line according to a Poisson law. However we verified 
that a regular spacing distribution does not modify the 
long time dynamics of the system. 

The time dependent velocity distribution has been 
monitored yielding the following results. First, we plot- 
ted on a double logarithmic scale the number of particles 
(respectively the root mean square velocity of a particle) 
as a function of time. The exponents a (respectively /3) 
are then extracted from the slopes of the curves. Lin- 
ear regressions were made for various sets of points. The 
retained value is the one corresponding to the best fit 
and the error is given by the maxinum deviation. Typ- 
ically, the fits were made for time intervals such that 
10^"^ ^ t 1. For longer times, the system starts to feel 
the boundary conditions. The lower limit of the time in- 
terval corresponds approximately to the beginning of the 
linear regime. Remark however that at this time about 
98% of the particles have already disapeared. The results 
are given in Table || and Figures ||-^. Data were averaged 
over 10^ samples for the Gaussian distribution and 10'^ 
for the uniform and Lorentzian cases, respectively. 




FIG. 2. Plot of the relative density (circle) and the root 
mean square velocity (plus) as a function of time, for a Gaus- 
sian initial velocity distribution ^(w) = exp(— w^/2)/%/27r. 
The periodic chain of length 2 initially contained 2^* par- 
ticles. Data are averaged over 10** samples. The two straight 
lines are obtained by linear regression over a subset of points 
(typically for t between 10~^ and 1). 



TABLE L Numerical values of the exponents a and (3 ob- 
tained from the simulations (see Fig. |^-^). The error is ±0.005 
for each value of a and /3. Results are given for three different 
initial velocity distributions. 



Velocity distribution 


a 


13 


a + 13 


Gaussian 


0.785 


0.195 


0.980 


Uniform 


0.795 


0.195 


0.990 


Lorentzian 


0.780 


0.195 


0.975 




FIG. 3. Plot of the relative density (circle) and the root 
mean square velocity (plus) as a function of time, for a uni- 
form initial velocity distribution 4>(v) = S(\/3— \v\)/\/l2, so 
that Af(0) = 1. The periodic chain of length 2 initially con- 
tained 2^^ particles. Data are averaged over 10"^ samples. 
The two straight lines are obtained by linear regression over 
a subset of points (typically for t between 10"'' and 1). 




FIG. 4. Same as in Fig. |3| for a Lorentzian initial velocity 
distribution (f>{v) = n~^{l + v^) 
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We see that the scahng law a + /3 = 1 is weU satisfied 
for the different distributions. For the Lorentzian one 
however, the results are less precise. This is due to the 
fact that the scaling regime sets in for longer times, and 
longer simulations would be needed to reach the same 
precision. 

Note that for a imiform velocity distribution, Ben- 
Naim et al. g obtained a = 0.76 and (3 = 0.22, but 
with a slightly different algorithm (which involved dif- 
fusion in addition to the ballistic motion) and a poorer 
statistics. 

From the values quoted in Table || we immediately 
deduce, at least for the three different velocity distribu- 
tions considered here, that the exponents a and (3 are the 



same. In view of the arguments presented in Section [I] 



we are led to conjecture that there is universality for a 
wide class of initial conditions. However, this universal- 
ity should not hold for initial distributions which diverge 
or vanish at w = [see Eq. (pO|)]. 

The scaling functions ^!{u) have also been measured 
for the three different initial distributions. In Figs |^-|^ 
the scaling functions ^'(u) are plotted versus the reduced 

— 1/2 

velocity u = v[M{t)] . For the Gaussian case, one sees 
(Fig. ^ that 5* keeps its Gaussian shape until t < 0.1. 
In fact a finer analysis of the six first even moments of 
* ((w^") with n = 0, 1, . . . , 5) shows that 'I' looses its 
Gaussian character when t > 0.2, i.e. when less than 50 
particles remain in the system. Ultimately, when less 
than 10 particles remain in the system, ^ tends to a bi- 
modal distribution. We emphasize again that this late 
stage behavior is an artifact of the finiteness of the sys- 
tem. In the thermodynamic limit, the true asymptotic 
behavior would be Gaussian like. Similar conclusions can 
be drawn for the uniform and Lorentzian cases (see Figs 
^ and 1^). Indeed, after a transient regime, ^' adopts a 
Gaussian profile until t > 0.2 . These conclusions are 



0.08 



0.06 



0.04 



0.02 



0.00 



: f=10 
□ f=10 
of=10 

f=10 
+ f=lO 

f=10 
a=1 



-4.0 -2.0 



0.0 
--vM 



2.0 



4.0 



U=VM '"(f) 

5. Plot of the scaling function 'i' for a Gaussian initial 



FIG. 



velocity distribution (see Fig. ^) as a fonction of the reduced 
velocity u = v[M{t)]~^^^ for 7 different times. Data are aver- 
aged over 10'' samples. One remark that *!/ keeps its Gaussian 
character as long as t < 1. 



again confirmed by the analysis of the moments. Thus, 
for these three distributions, we conclude that beside the 
universal behavior of the exponent, the scaling functions 
^{u) are also the same. There is an attractive Gaussian 
like scaling distribution in the long time regime: 



*(m) 



/27r 



exp(-uV2)- 



(35) 



Here again, we are led to conjecture that such a behavior 
would be valid for any velocity distribution which takes 
a finite non-zero value ai v — and which is regular near 
w = 0. 

In addition, the range of validity of our scaling postu- 
lates can be tested. Indeed from Eqs (17p^), one sees 
that 



0.08 - 



0.06 



0.04 



0.02 



0.00 



f * 



ot=^o 

+ f=10"' 



-4.0 -2.0 



2.0 



4.0 



0.0 

M=vM"'(f) 

FIG. 6. Plot of the scaling function >I/ for a uniform initial 
velocity distribution (see Fig. ^) as a fonction of the reduced 
velocity u = v[M{t)]~^^^ for 5 different times. Data are av- 
eraged over 10"^ samples. One remarks that after a transient 
regime not represented here, 'i' is attracted towards a Gaus- 
sian like scaling distribution (as soon as f > 10~*). 



0.08 - 



0.06 - 



0.04 - 



0.02 - 



0.00 



,f=10" 
f=10"' 

<t=^o'' 
^t=^ 



4 



4 



-4.0 



2.0 



4.0 



-2.0 0.0 

M=vM"'(f) 

FIG. 7. Same as in Fig. |^ for a Lorentzian initial velocity 
distribution. 
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n{t) = const (n{t)[M{t)]^^'^y, (36) 

1 /2 

SO that the log-log plot of ri(i) versus n(t)[Af(i)] ' should 
give a straight line. In Fig. ^ we reproduce this plot for 
a Gaussian initial distribution. We remark that Eq. ( |3^ ) 
is satisfied for times as short as i ~ 10~^ (at this time 
less that 20% of the particles have already reacted) and 
as long as t ~ 5. Note that the slope (i.e. the exponent 
a) obtained by linear regression is slightly larger than 
the value obtained from Fig. ^. This comes from the 
fact that the sum a + (3 given by the value of the expo- 
nents measured on Fig. ^ gives 0.980 rather than strictly 
1. Multiplying the slope of Fig. js] by 0.980 reproduces 
indeed the value of a quoted in Table |. For the uni- 
form and Lorentzian distributions the range of validity 
of the scaling postulates is smaller, beginning only near 
t ~ 10~^. This fact confirms the particular role played 
by the Gaussian distribution in this problem. 




FIG. 8. Log-log plot of the density as a function of 
n{t)[M{t)]^^^ for a Gaussian initial velocity distribution 
(data are the same than Fig. ^, note that large values of 
n{t)[M{t)]^^^ correspond to short time). The comparison 
with a straight line is excellent both for short and long times. 
This shows that the range of validity of our scaling postulate 
is very wide for a Gaussian initial condition. 

V. CONCLUDING REMARKS 

We have studied the kinetics of ballistic annihilation 
for a one-dimensional ideal gas with continuous velocity 
distribution. Starting from an exact analytical approach 
previously derived, we established a scaling theory for the 
long time behavior of such systems. The validity of this 
scaling theory has been tested numerically for three dif- 
ferent initial continuous velocity distributions (Gaussian, 
uniform and Lorentzian) . Both the dynamical exponents 
and the scaling functions are the same for the three cases. 
This led us to conjecture that all the continuous velocity 
distributions which take a finite nonzero value aA, v — 
and which are regular near v = are attracted towards a 
Gaussian distribution and thus belong to the same uni- 
versality class. Despite several attempts we have not yet 



been able to prove this conjecture in the framework of 
our exact analytical approach. We have only shown that 
a Gaussian distribution is compatible, in the long time 
regime, with the exact dynamical equation (^). How- 
ever, we have not been able to prove that this Gaussian 
distribution was the only possible solution. 
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APPENDIX: NUMERICAL ALGORITHM 

To simulate the ballistic annihilation in one dimension, 
the simplest algorithm is probably the standard molec- 
ular dynamics: starting from a given configuration, one 
identifies the shortest collision time, removes the two col- 
liding particles and calculates the new positions of the 
remaining particles at this time. Starting from this new 
configuration, the process is iterated. This algorithm is 
very simple but not very efficient, the computing time 
increasing with the number of particles N as 0{N'^). 

The numerical algorithm we used instead has been 
largely inspired by the one developed by Krapivsky et 
al. |4j . The idea is to establish the list of all the ( "true" ) 
collision times arranged in chronological order. From the 
initial condition we compute the collision times of each 
particle with its right nearest neighbor and sort those 
times into an ascending series, called A, using a stan- 
dard sorting algorithm (see for example |jl^). The short- 
est time of this set corresponds to the first "true" colli- 
sion time and is the first member of the list of the "true 
collision times" ; simultaneously it is removed from A. 
Then one removes the pair of particles (n,n + 1) corre- 
sponding to this first collision. As a consequence, the 
collision times associated with the pairs (n — 1 , n) and 
(n + l,n + 2) should be discarded from A, producing a 
truncated sorted list called A' . The collision time of the 
new nearest neighbor pair (ti — 1, rt -I- 2) is computed and 
is the first element of a new unsortcd list, called Af. The 
process is then iterated starting with the sorted list A' 
as long as its first element is shorter than the shortest el- 
ement of J\f. When this is no longer true, we merge both 
lists {A' and J\f) into a new one which is sorted. This 
last list replace A and the process continues until at most 
one particle remains in the system. For the continuous 
velocity distribution we considered, this merging step be- 
tween A' and M occurs very rarely. For example, for a 
Gaussian initial velocity distribution, it takes place ap- 
proximately hundred times for 2^^ iterations of a system 
containing initially 2^® particles, and the whole simula- 
tion used about thirty seconds of CPU on a HP 9000 Serie 
700 workstation. The computing time increases with the 
number of particles N roug hly as C'(Ar5/4 ^^j^y 
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